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ABSTRACT 

Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass 
can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds 
(LMC and SMC). However, it would be straightforward to interpret such an estimate only if the 
properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. 
In this study we explore whether the environment of the Milky Way affects the properties of the SMC 
and LMC such as their velocities. To test for the effect of the environment, we compare velocity 
distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and 
in the second sample of such halos selected with additional restrictions on the distance to the nearest 
cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). 
We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as 
compared to the halos of similar mass in the sample without environmental constraints. For example, 
the fraction of subhalos matching the velocity of the LMC is 23±2% larger in the LG-like environments. 
We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in 
the control sample and find that the environment does not significantly affect the derived likelihood. 
We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass 
log(M2oo/M0) = 12. 061° 19 (90% confidence interval). We also explore the incidence of close pairs 
with relative velocities and separations similar to those of the LMC and SMC and find that such pairs 
are quite rare among ACDM halos. Only 2% of halos in the MW mass range have a relatively close 
pair (Ar < 40kpc and As < 160 km s~ x ) of subhalos with circular velocities w C i rc > 50 km s" 1 . Pairs 
with masses and separations similar to those of the LMC and SMC (AfMC = 23.4 ± 10 kpc, and 
Asmc = 128 ± 32 km s^ 1 ) are found only in one out of 30000 MW-sized halos. Interestingly, the 
halo mass likelihood distribution for host halos constrained to have MC-like close pairs of subhalos is 
quite different from the global likelihood from which the MW halo mass constraint discussed above was 
derived. Taking into account the close separation of the MCs in the Busha et al. 2011 method results 
in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately 
by a factor of two. 

Subject headings: Galaxy: fundamental parameters, halo — galaxies: Magellanic Clouds — dark 
matter 



1. INTRODUCTION 

In the CDM scenario the Milky Way (MW) halo 
formed by accretion and disruption of smaller ha- 
los, and some of them survived this process as self- 
bound substruc tures or subhalos orbiting around the 
hos t halo (e.g.. IKlypin et al.l I1999U IMoore et all 119991 
see iKravtsovl 120101 for a recent review). The major- 
ity of these subhalos are likely devoid of stars, but 
some of t hem are massive enough to host dwar f galax- 
ies (e.g.. IKauffmann et all Il993at IKlypin et al.l Il999bt 
Bullock et al.l 120001: IKravtsov et all I2004t iFont et ail 
20111) . The MW and its satellites offer a unique labo- 
ratory for testing CDM predictions for halo formation 
and associated substructure, because properties of many 
of the satellites have been studied in details in observa- 
tions, including proper motions and rotation curves (e.g., 
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IKallivavalil et~alll2006blial) . 

It is still debated how typical is the satellite popula- 
tion of the MW, but both theoretical models and observa- 
tions indicate that the incidence of the Magellanic Clouds 
(MCs) satellites is rare for MW-sized galaxies. Semi- 
analytic galaxy formation models predict that only ~ 
10% of the MW -sized galaxies have satell i tes as bright as 
the LMC (e.g.. [Kauffmann et al.| | 1993bt iKoposov et al.1 
I2009fl . iBovlan-Kolchin et al.l (|2011h used large sample of 
MW-sized halos extracted from the Millenium-II simula- 
tion and found that 20% MW-mass halos host an LMC 
or SMC, and only ~ 2.5% host bq t h MCs . A similar fre- 
quency was found bv I Busha et al.1 ()201 lbf ) using the Bol- 
shoi simulation of the concordance cosmology and abun- 
dance matching ansatz to assign stellar masses to ha - 
los and subhalos. In observations. Uames fc Ivory) ([201 ID 
searched for star-forming satellites around 143 luminous 
spiral galaxies and found that two-thirds of central galax- 
ies have no satellites down to luminosity and star forma- 
tion rates well below those of the MCs . Using the Sloan 
Digital Sky Survey (SDSS) , EhleLalJ §20uJ) computed 
the occurrence of the satellites similar to the MCs in lu- 
minosity around MW-sized galaxies and found that 11% 
have one, and only 3.5% have two MCs within a radius 
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of 150kpc. Similar values for the freque ncy of MCs are 
found in the SDSS in a recent study by iTollerud et al 



2011 ) , and in the GAMA survey by iRobotham et al 



20121 ) . What makes Magellanic Clouds even more pecu- 



liar is the fact that they are likely to be an interacting 
pair (see § l4.3l below) and have rather high velocities with 
respect to the Milky Way. 

One caveat in this debate is that some satellite prop- 
erties depend sen sitively on the MW halo mass (e.g., 
iWang et al.ll2012T ). Although there are sev eral methods 
to constraint the MW hal o mass (see e.g., iBusha et al.l 
1201 lat IDeason et al.l 120121 and references therein) , it is 
still quite uncertain. In this study we focus on the 
methods that use satellite properties to c onstrain the 
MW ha lo mass. In particular, as shown by IBusha et al.l 
(|2011al Bll hereafter), the MW host halo mass likeli- 
hood can be computed using the observed MCs prop- 
erties and a statistically representative sample of CDM 
halos to evaluate the likelihood that a given halo would 
have each or all MCs properties, such as circular velocity, 
host halo distance, and velocity within the host. Such ap- 
proach provides a new interesting constraint on the MW 
virial mass, but there are certain questions that need to 
be addressed to assess its reliability and interpretation. 
One of such questions is whether the environment of the 
Local Group biases properties of the MW satellite popu- 
lation, and MCs in particular. We explore this and other 
related questions in more detail in this paper. 

Another important question is the origin and dynam- 
ics of the MCs. The accretion history of the MCs is 
still not clear, the presence of the Magellanic Stream 
(MS hereafter), a filament of gas extending 150° across 
the sky, with an apparent spatial and chemical associ- 
ation with the MC s, is interpreted as a tidal tail (see 
IBesla et al.l I2012al and references therein) . There are 
several clues which indicate that t he MCs are bound , 
and have been inte racting recently ([Diaz & Bckki 2012; 
IBesla et all l2012bT) . and there is also some evidence of 
the SMC stars accreted onto t he LMC (|Qlsen et al.l 
2011; Cas etti-Dinescu et al.ll2012|) . Cl ose satellite pairs 
are r a re in nearby MW-size d galaxies (jRobotham et al] 
I2012t IJames fc Ivorvl 1201 ID. Proper motion measure- 
ments ioTti^~MCs<^^^^^MMl^DQMi) indicate 
that they have high velocities not aligned with the Mag- 
ellanic Stream, and suggest that the MCs could be 
bound and on their first or second pericenter passage 
into the MW (IBesla et al.l [20071 l2010l 1201 2bD . Other 
studies also sugg est that the MCs were accreted in 
the s ame system (|D'Onghia fc Lake! 120081 iNichols et all 
1201 If) . Bll estimated from the Bolshoi simulation that, 
for a typical MW-sized host, there is a ~ 72% proba- 
bility that the MCs were ac creted within the l a st Gyr 
IBovlan-Kolchin et al.l ([201 ID and IBesla et al.l (|2012aD 
also favor the late accretion scenario for the MCs. 

The paper is organized as follows. We describe the 
simulation and halo catalogs used in our study in § [21 
while in §[3] we describe the data samples for our MW-, 
MCs-, and LG- analogues. We presents the MW halo 
mass likelihood in different environments and under dif- 
ferent assumptions on whether MCs are independent ve- 
locity samples in § @] and § 14.31 We present discussion 
and our conclusions in § [SJ 

2. SIMULATION AND HALO CATALOGUES 



To carry out our anal ysis we use halos ex tracted from 
the Bolshoi simulation ([Klypin et al.l 120111) , which fol- 
lowed evolution of 2048 3 particles in 250/i _1 Mpc cu- 
bic volume assuming concordance flat ACDM cosmology 
with parameters: J7 m = 1 — — 0.27, ilb = 0.0469, 
h = iio/(100) = 0.7, o 8 = 0.82 and n s = 0.95, compati- 
ble with combined cons traints from WMAP, BAO, SNe, 
and cluster abundance ([Komatsu et al.ll2011h . The high 
spatial and mass resolution and relatively large volume 
make the Bolshoi simulation well suited for providing a 
base halo sample for our study. 

The halo catalog we used was c onstructed using the 
Bound Density Maxima h alo finder ([Klvpin fc Holtzmanl 
ll997tlKlvpin et al.l[l999i) . which identified 2285977 ha- 
los down to the resolution completeness limit of w C i rc « 50 
km s _1 (corresponding to « 110 particles or M200 ~ 1-7 x 
10 10 Mq). The circular velocity 

^circ is the maximum 
of the circular velocity profile ?; c i r c(r) = ^JGM(< r)/r. 
The parameters of the halo finder were set such that 
the density maxima are not allowed to be closer than 
10/i _1 kpc, and the finder keeps only the most massive 
density maximum if that happens. Halo center is identi- 
fied with the particle location which has the largest local 
density, and bulk halo velocity is computed as the av- 
erage velocity of the 30 closest neighbors of the central 
particle. The algorithm computes a number of halo prop- 
erties after an iterative procedure to remove unbound 
particles. 

Throughout this paper we defined halo mass, M200, as 
the mass enclosed in a sphere of radius -R200 with the 
density 200 times the critical density of the universe at 
the rcdshift of analysis. Another common mass definition 
is M v i r , the mass within the radius enclosing the mean 
overdensity of 358 with respect to the mean density of the 
universe (or overdensity of 358 x 0.27 ~ 97 with respect 
to the critical density) ([Bryan fc Nor man 1998). In the 
Bolshoi simulation, we find a relation of M v - lr / 'M200 = 
1.21 for host halos of v C i IC « 220 km s . 

We use MW-sized host halos in the Bolshoi catalogs 
to search for subhalos with velocities and positions sim- 
ilar to those of the MCs. However, positions of satellites 
change on short timescales due to their motion along 
their orbits. Therefore, to increase statistics, we stack 
the halo catalogs of several simulation snapshots close to 
z = separated by Aa = 0.003 (~ 42Myrs at z = 0). 
During the time interval between snapshots a typical 
MC should move ~ 15kpc along its trajectory. We 
stack satellites of the last 30 snapshots, and the total 
time difference between the first and last snapshot is 
~ 1.3Gyrs (or Az < 0.1), so we can neglect any evo- 
lution effect in MW mass halos (iMunoz-Cuartas et al.l 
[Mil iDiemer et al.ll2012T: ICuesta et al.H2008D . 

For a given MW-sized halo at z = we should, in 
principle, consider the last N simulation outputs to trace 
trajectories of all satellites and check if they match MCs 
contraints at some time in the recent past. However, this 
is computationally expensive. Instead, in our analysis we 
consider each snapshot as an independent realization of 
halo properties, which effectively increases the simula- 
tion sample by a factor of N. In this approximation we 
neglect any correlation in the positions and velocities of 
satellites between snapshots. For example, in the ex- 
treme case of a purely circular orbit, a single LMC or 
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SMC analogue would be counted N times. However, for 
realistic eccentric orbits such double counting is quite 
rare, especially for satellites with relatively small radial 
distances to host centers similar to those of the SMC 
and LMC. On the other hand, we have a large number of 
hosts and satellites in each snapshot, and the randomness 
of their orbital configuration produces consistent distri- 
butions and average fraction of satellites matching any 
given set of constraints, for any number of snapshot se- 
lected. We have tested for the double counting effects us- 
ing different number of snapshots, N, computing the dis- 
tribution of satellites matching several set of constraints, 
and we have found no significant differences in the prop- 
erties of satellites for TV < 30. 

3. HALO SAMPLES AND OBSERVATIONAL CONSTRAINTS 

In our analysis we use the following three main obser- 
vational measurements. 

3.1. The Milky Way halo mass 

To select MW analogues from the Bolshoi-derived halo 
catalogues, we select halos within the broad mass range 
M 2 ooc = 0.8 — 2.9 x 10 12 M Q , which covers the range 
of current ob servational constrain ts: e.g., using HI gas 
distribution dKalberla et al.l 120071). kinematics of stars 
dSofue et all 120091: [Xue et all 120081: iGnedin et all |2010t 
iWilkinson fc Evansl 119991: iDeason et al.1 12012ft . satel- 
lite dynamics dWatkins et al.1 l2010h . esc ape velocity 
(iSmith et ail 120071 ). and timing argument ()Li fc White! 
120081 ). There are (~ 57000) host halos in this mass 
range in the Bolshoi simulation at z = 0, which contain 
~ 115000 subhalos with t> C i rc > 50 km s . 

3.2. Magellanic Cloud analogues 

We follow iBusha et al.l (|2011aD and use the following- 
observed properties of the Magellanic Clouds to select ap- 
propriate MC analogues among subhalos in the MW halo 
analogues and to constrain the halo mass of the Milky 
Way: the distance to the host center ro, the total speed 
relative to the host center s, and the subhalo circular ve- 
locity fcirc- F or these quantities we us e recent HST mea- 
surements bv lKallivavalil et al.l (f2013h : 76.1 ±7.6 
kms" 1 , r = 50±5kpc, s = 321±24 km s~Hot the LMC, 
and w c i rc = 60±5 km s _1 , ro = 60±5kpc, s = 217±26 km 
s _1 for the SMC. The halos hosting the MCs analogues 
are required to have at least two subhalos with v C j rc > 50 
km s _1 and we will consider only the two subhalos with 
the largest circular velocities^. 

Note that ro errors are inflated from their actual ob- 
servational values to improve statistics of the sample of 
the MC analogues. The range of ro we use corresponds 
to typical radial displacement of subhalos along its orbit 
between consecutive snapshots. These parameters are 
used everywhere in this paper, except in § 4.2, where we 
use a different definition for the MCs analogues assum- 
ing a single bound system, and in § 5, where we include 
two additional constraints: the relative separation and 
velocity of the clouds. 

5 In Bll, the MCs analogues are selected in hosts that have 
exactly two subhalos with v c - lIC > 50 km s , we will explore the 
effects of this difference in the Results section. 



3.3. The Local Group analogues 

The Milky Way is not an isolated g alaxy, but is located 
in a pair with M31 (~ 7 70 kpc awav iMcConnachie et al.l 
l2005t iRibas et"aL1 120051 ) and is surrounded by a num- 
ber of smaller galaxies, collectively known as the Local 
Group of galaxies. On larger scales the environment of 
the Local Group is rather low-density: in a sphere of 
50 Mpc radius around the L G, the estimated de nsity is 
~ 3 times lower than average (|Karachenstevll2012l ) , while 
in a sphere of 5 Mpc the density around the Local Group 
is close to the mean densi ty of the universe (|Klypin et al.l 
120031: lKarachentsevll2005t ) . In addition, the nearest clus- 
ter to the Local Grou p is the Virgo Cluster ~ 16.5 Mpc 
away (|Mei et al.ll2007D . 

To explore whether environment of the Local Group on 
different scales affects statistics of the MC analogues, we 
derive several MW analogue samples that mimic differ- 
ent aspects of the real MW environment: the host halos 
in the P sample have an M31-sized companion in rel- 
ative isolation with no other large neighbor; host halos 
in the LGP sample are a subset of halo pairs from the 
P sample, but with additional environmental constraints 
designed to more closely mimic the LG environment; fi- 
nally, host halos in the S sample include all host halos 
that are not in included in the LGP- and P- samples. 
These sample definitions and naming convention will be 
used for both host halos in the MW halo mass range 
(Section 3.3) and halo samples in a wide range of masses 
(Section 4). Note that for the P and LGP samples the 
mass of the M31-like companion is always fixed to the 
same mass range, even as its MW analogue halo mass is 
varied within a wider range. 

Below the sample definitions are described in more de- 
tail. 

The P sample: the host halos are required to have 
a companion similar to M31, which we define using the 
same mass range of Mima = 0-8 — 2 .9 x 10 12 M Q (e.g., 
iWatkins et"aB2010tlEvans et al.ll2000D . The pairs are re- 
quired to be relatively isolated and not be part of a triplet 
or a larger group. As a quantitative isolation criterion 
for the pair we use the force constraint F i]Com < nFi 2 , 
where Fi iCorn is the gravitational force between the pair 
and any neighbor halo i within a 5/i~ 1 Mpc radius of the 
pair center-of-mass, and Fyi is the force between halos 
in the pair, and k is a constant parameter. The isolation 
criterion becomes increasingly strict for decreasing val- 
ues of k. The Milky Way and M31 do not have massive 
neighbors within 5 Mpc, and should thus have n < 0.1. 
The actual value of k is, however, uncertain, and we ex- 
plore a range of values in the further analysis. We use 
K = 0.1, 0.15, 0.2, 0.25 which results in 205, 378, 598, 
810 LG-like halo pairs, respectively. From the a poste- 
riori analysis of satellite properties, results are largely 
consistent for different k values and we thus use k = 0.2 
in this paper to maximize the statistics of the MW hosts. 
The 598 LG-like halo pairs we found using n = 0.2 repre- 
sent ~ 2% of all halos in the MW mass range, and ~ 80% 
of these pairs are bound under the two-body point mass 
energy approximation. 

The LGP sample of halos is a subset of halos in 
the P sample with additional constraints on the local 
and global environment to more closely mimic the envi- 
ronment of the Local Group. Namely, we require that 
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Fig. 1. — Distribution of the local galaxy overdensity (left panel), distance to the nearest cluster-sized (M200 > 1-5 X 10 14 Mg) halo 
(center panel), and distance from the pair companion for halos in the P sample (solid black line), LGP sample (short-dashed line), and 
S-sample (red long-dashed line). The maximum local overdensity, the minimum distance to the closest cluster, and the minimum pair 
separation used as additional constraints for the LGP sample of the MW analogues in the LG-like environments are shown by the black 
arrows. The MW is ~ 16.5Mpc away from the Virgo cluster (middle panel, red arrow) and separation between the MW and M31 is 
~ 770kpc(right panel, red arrow). 

are w 1.3% of the total number of host halos in the MW 
halo mass range. This indicates that the LG environment 
of the Milky Way is rather rare for halos of this mass, the 
fact also indicated by the "coldness" of t he local velocity 
field of galaxies (e.g.. iKlvpin et al1l2003[ ). 

The figure shows that the P sample halos are located 
in environments with similar overdensity distribution to 
the S-sample halos, but have a narrower distribution due 
to the force constraint and the chosen k = 0.2 value, 
which eliminates responsible pairs close to other larger 
halos and single isolated halos. 



halos in this sample do not have a cluster-sized halo 
more massive than M200 = 1.5 x 10 14 M Q within 12 Mpc. 
This mas s limit is close to the mass of the Virgo Cluster 
M (e.g., iFouaue et all [20011 iNulsen fc BohringeH H995l 
and references therein). For the local environment we 
compute the galaxy density field using Voronoi tessella- 
tion (VT hereafter) o n halo positions and masses (sim- 
ilar to the method of [Schaap fc van de Weveaert| [2000; 
iGonzalez fc Padillal [2010f>. The VT partitions the vol- 
ume into cells, where each cell is associated with a single 
host halo. The shape and volume of each cell is defined by 
distribution of halo neighbors. The adaptive local den- 
sity can be computed using the local cell volume around 
each halo and the enclosed halo mass, but instead we use 
also neighbors cell volumes and masses to compute the 
average, in this way we generate a smooth density field 
where the typical number of direct neighbors around a 
halo is w 14. We use only the host halos with mass 
higher than M 2 oo = 1.5 x 10 10 M for density computa- 
tion, and define the local overdensity as 8 = (p— (p))/{p), 
where (p) is the mean of the density distribution. The 
local density of the LG is not well constrained, but the 
abundance of luminous galaxies within 3 — 6 Mpc is close 
to the average densi ty of galaxies in the local universe 
(jKarachentsevi l2~005h . Another local environment con- 
straint we include is the distance to the M31-like pair 
companion. 

Figure [T] shows the distribution of the local density, dis- 
tance to the closest Virgo-sized halo, and distance to the 
M31-sized pair companion for the MW analogue halos in 
the P sample, and the cuts we impose to define the LGP 
sample. Specifically, we select only halos with the local 
overdensity smaller 5 = 1.04 to avoid systems located in 
the highest density decile. The mean overdensity for the 
S-sample is located at log (5 + 1) ~ 0, and for the LGP 
sample it is located 0.5 dex lower. 

We exclude halos with Virgo-sized neighbors closer 
than 12 Mpc and require that distance to the pair com- 
panion is in the range 0.4 < Ar < 1.2Mpc to avoid close, 
possibly merging pairs, but to include pairs with sepa- 
rations similar to the actual distance between the Milky 
Way and M31. These constraints eliminate about a third 
of halos from the P sample; the halos in the LGP sample 



4. THE MW HALO MASS ESTIMATE 

In the context of the Bll method, in order to estimate 
the halo mass of the Milky Way from the properties of 
the MCs, such as circular velocity, velocity and position 
relative to the center, we explore two aspects which can 
affect the mass estimate: 1) the environment, in particu- 
lar whether differences in the environment correspond to 
the differences in the subhalo populations of halos, and 
2) whether it matters if MCs are treated as two indepen- 
dent dynamical samples or a single tracer (a bound pair 
sharing a common translational motion of their center of 
mass) . 

4.1. Effect of the environment 

To test for possible effects of the MW environment on 
subhalo statistics, we extract all subhalos with v c - n - c > 50 
km s _1 , and compute their distance to the host center 
ro and the total speed relative to host s. Subhalos very 
close to the host center (ro < 20/i _1 kpc) are removed to 
avoid artefacts or resolution-related problems in the halo 
identification procedure. 

Figure shows abundance of subhalos with w circ > 
50, 60, 70 km s in four mass ranges around 10 12 M Q . As 
expect, the average number of subhalos increases with 
increasing host halo mass and decreases for increasing 
subhalo circular velocities. The figure shows that there 
is no significant difference in the abundances of subhalos 
in hosts of the S and LGP samples, which means that 
the Local Group environment does not appreciably af- 
fect abundance of the massive subhalos in the MW-sized 
hosts. 
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Fig. 2. — Average number of satellites in the S and LGP samples 
for halos in different mass ranges around most likely range of the 
MW halo mass, and satellite circular velocity ranges defined in 
legend. The error bars of the mean values are computed using the 
jackknife method. Top-right legend quotes the fractions of halos 
with Ng\i~ members for subhalos with v c i IC > 50 km s" 1 , and host 
halo masses in the range 0.8 < A/200 < 2.9 X 10 12 M Q . 



0.004 



> 0.003 



0.002 



0.001 



0.000 



S-somple 
LGP-sample 




100 



200 300 400 
s [km s _, l 



500 



Fig. 3. — Distributions of the total subhalo velocity in the frame 
of their host halo for hosts in the S and LGP samples. Blue arrows 
indicate the velocities of the SMC and LMC. The uncertainty of the 
distribution for the LGP sample is shown by the dashed lines (es- 
timated using the jackknife method). The environment has small, 
but statistically significant effect on the distribution of velocities: 
the fraction of subhalos with the LMC velocity in the LGP sample 
is 23 ± 2% higher than in the S sample. 

Figure [3] shows the distributions of the total subhalo 
velocity in the frame of their hosts for halos in the S 
and LGP samples. The fraction of satellites with veloc- 
ities comparable to that of the LMC is 23 ± 2% higher 
in LGP sample than in the S sample (there is no sig- 



nificant change in the distribution at the SMC speed). 
Therefore, the fraction of satellites matching the veloc- 
ity of the LMC is somewhat enhanced in the LG-like 
environments. This can affect the MW mass estimate 
because LMC velocity is the main constraint driving the 
mass likelihood distribution to larger masses (see Bll 
and Fig. [5] below). 

We have tested whether the choice of k parameter used 
in the isolation criterion (see § !3-3[) influences the veloc- 
ity distribution and have found that the result velocity 
distribution is almost the same for values of k in the 
range 0.1 < n < 0.25, with a weak trend towards higher 
fraction of subhalos in the LGP sample matching the ve- 
locities of the LMC and SMC for lower k values. The 
magnitude of the trend, however, is comparable to the 
uncertainties in the velocity distribution. 
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Fig. 4. — Fraction of subhalos matching the velocity of the SMC 
and LMC in the S (black symbols and lines) and LGP (red sym- 
bols and lines) samples for three different halo mass ranges around 
M200 = 1O 12 M0 (the first three connected bins) and the total mass 
range used to define the MW halo analogues (the gray shaded re- 
gion). 

In figure 2] we show the fraction of subhalos with total 
velocities relative to the host center similar to those of 
the LMC and SMC for host halos in different mass ranges 
around M200 = 10 12 M Q . There is a small enhancement 
in the frequency at the LMC velocity, but no difference 
at the SMC velocity. In the LGP sample the fraction of 
host halos with subhalos matching the LMC velocity is 
larger by 29 ± 10%, 16 ± 4%, and 10 ± 2% than in the S 
sample for 0.8 - 1.2, 1.2 - 1.9, and 1.9 - 2.9 x 10 12 M Q 
mass ranges respectively. 

We present the likelihood distribution for the MW halo 
mass from the constraints from the Magellanic Cloud 
properties in Figure [5] To compute the likelihood we 
used the S and LGP sample constructed without any 
restriction on the M200 mass of halos. In addition to the 
likelihood for all of the constraints combined for S-sample 
only, we also show likelihoods resulting from using only 
one of the constraints for both S and LGP samples. In 
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Fig. 5. — Likelihood distribution for the MW halo mass (black 
line), M20O1 m the S sample based on the properties of the LMC 
and SMC allowing the full halo mass range instead of the range 
A/200 = 0.8 - 2.9 X lO 12 Af . Lines of different colors show likeli- 
hoods when only one constraint is used: exactly 2 subhalos with 
v c > 50km/s (blue line), circular velocities of the MCs (red line), 
radial positions of the LMC and SMC (orange line), total velocity 
relative to the host center of the LMC and SMC (green line). Dot- 
ted lines are corresponding distributions for the LCP sample using 
the same color key, but we omit the distribution using all con- 
straints together in this case due lack of matching systems. The 
MW mass estimate in the S sample is log(M 2 oo/M ) = 12.06to'ig 
(90% confidence interval). 

particular, we present the likelihood for the constraint 
that host halos have exactly two subhalos with v c - n - c > 50 
km s with all other constraints (circular velocity, radial 
distance, and total velocity) applied after this condition 
to compare to the results of Bll. 

We chose not to enforce the two satellite condition in 
our estimate of the likelihood, instead we allow the host 
halo to have any number of satellites to increase the sam- 
ple size, but we will consider only the two largest ones. 
The mass estimates obtained with and without enforcing 
the two satellite condition are presented in Table [1] The 
final mass estimate using all of the constraints together 
is only computed for the S sample, in which there are 
40 MC-like satellite pairs. There are only 3 such satel- 
lite pairs in the LGP sample, and therefore we do not 
attempt to derive the total likelihood for this sample. 
However, given that the likelihood distributions for indi- 
vidual property constraints are similar for the S and LGP 
samples, we expect that the total likelihood is similar for 
the combined constraints as well. 

The MW halo mass estimate for the S sample is 
log(M 2uO /M ) = 12.06+^;^ (90% confidence interval), 
which is in general agreement with the Bll result in the 
central value. We present results for different samples 
and constraints in Table [IJ where we also include both 
the 68% and 90% confidence interval errors. Due to sig- 
nificant deviations of the likelihood from the log-normal 
form in the tails, the 90% errors are significantly larger 
than the 68% ones and we therefore chose to quote the 
90% errors. 



Note that there is a key methodological difference in 
the way the likelihood distribution was evaluated in Bll 
and in our analysis. In our mass estimate we use all the 
constraints together simply as the distribution of halo 
properties that satisfy the constraints. This way any 
correlations between properties, such as expected corre- 
lation between radial distance to the host center and total 
velocity are taken into account automatically. Bll, on 
the other hand, assume that the probability distributions 
for each constraint are independent (see their eq. 3). If 
we multiply the probability distributions for radial, speed 
and circular velocity constraints we get an estimate of the 
MW mass which should be similar to the method of Bll. 
For the S sample this gives log(M 2 oo/M ) = 12.24±g;|| 

and for the LGP sample log(M 2O o/M ) = 12.44+g;^. 
Thus, treating constraints as independent results in a 
small overestimate of the mass. 

4.2. Dependence on the SMC and LMC constraints 

The second key assumption of the Bll analysis is that 
the properties of the two Magellanic Clouds are not cor- 
related. However, the MCs are likely a bound pair, and 
their velocities and radial positions can thus be expected 
to be correlated. We have tested whether this assump- 
tion affects the MW mass estimate by comparing results 
in the following two cases: 1) properties of the MCs are 
treated as independent and we compare mass likelihood 
using both MCs, only the LMC, or only the SMC; 2) the 
MCs are considered to be a bound pair and we compute 
the likelihood using average properties of the pair, rather 
than properties of the two MCs independently. We ex- 
plore the latter scenario in the next subsection. 

For the case 1 we compute the mass likelihood for the 
LMC and SMC independently and find log(M 20 o/M Q ) = 
12.18lg' 2 3 usm g the LMC only, consistent with the result 
using both MCs discussed above. For the case when we 
use velocity of the SMC only, the mass is not well con- 
strained as the likelihood extends to considerably lower 
masses. Thus, the MW mass estimate is dominated by 
the properties of the LMC, as is expected since it is more 
massive and has a higher velocity. 

The fact that the main mass constraint comes from 
the satellite with the largest velocity is generic. If one 
considers distribution of absolute magnitude of satellite 
velocities as a function of the radial distance to the host 
center, at any given r the distribution is broad but has a 
sharp cut off at the velocity close to the escape velocity of 
the host. It is this sharp cut off that constrains the mass, 
and the constraint is due to the satellite with the largest 
absolute velocity at a given radius. This fact was recently 
used to constrai n the halo mass of the Milky Way with 
Leo I satellite by Bovlan-Kolc hin et al.l (|2012f K 

In principle, we can include additional properties of 
the LMC and SMC, such as as their separation or rela- 
tive velocity (~ 23kpc and 128 ± 32 km s _1 , respectively 
iKallivavalil et all l2006al 120131) , in the derivation of the 
mass constraint. However, this is difficult in practice be- 
cause very few host halos in the Bolshoi simulation match 
all of the properties of the LMC and SMC. We examine 
the incidence of the MC-like close pairs of satellites in 
the next section. 

Here we adopt a different approach, in which we as- 
sume that the LMC and SMC are a bound pair and 



Satellites in MW-like hosts 



7 



can be considered as a single subhalo. We then con- 
sider the pair as a single velocity tracer and use the 
velocity of the center-of-mass of the pair as a con- 
straint. We calculate the center-of-mass velocity as- 
suming that the MCs move in the same direction and 
neglecting peculiar velocity around the center-of-mass: 
||scm|| = + s 2 M 2 \\/(M 1 + Af a ) = 298 ± 52 km 

s _1 . We also treat the pair as a single subhalo with 



the circular velocity of v c 



The 



85.4±^ km 

masses are evaluated numerically from the M200 — ^circ 
relation for subhalos in the entire Bolshoi simulation. 
The asymmetry in the error ran ge is because we in clude 
larger values for the LMC v c j rc (Olsen ct al.ll20TTI ). Fi- 
nally, we adopt the radial distance of the center-of-mass 
of r cm = (nMi + r 2 M 2 )/(M 1 + M 2 ) = 54 ± 5 kpc. 
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Fig. 6. — Likelihood distribution of the Milky Way halo result- 
ing from the SMC and LMC being considered a bound system 
corresponding to a single subhalo (with parameters quoted in the 
text). The distribution for the S sample is shown by the black solid 
line, and for the LGP sample by the green solid line. The corre- 
sponding halo mass constraints are log(M2oo/MQ) = 12.13^q25 
and log(M 2 oo/M ) = 12.17±o'j4 (errors correspond to the 90% 
confidence level) for the S and LGP samples, respectively. Mass 
estimate for the S sample from Fig. 1 derived when properties of 
both MCs are used separately is shown for comparison by the dot- 
ted red line. 

Figure [6] shows the likelihood distribution for the MW 
halo mass estimated with such constraint for the S and 
LGP samples. We also include the mass estimate from 
figure [5] for comparison. The halo mass in this case 
is constrained to be log(M 20 o/M Q ) = 12.13±g;|| (90% 
confidence interval) in the S sample and log(A/2oo) = 
12.17l°;|° in thc LGP sample. We do not find any sig- 
nificant variation with environment in this case, and the 
mass estimate errors are consistent with the case when 
we constrain the mass using both MCs as independent 
tracers. The only difference is that the use of the MCs 
as a single tracer eliminates the tail of the likelihood to- 
wards small masses (M 2 qq ~ 10 11 M Q ), which is due to 
the fact that subhalos with circular velocities as high as 



u c irc = 85 km s _1 are highly unlikely in such small host 
halos. As noted above, the lower velocity SMC does not 
influence the constraint due to the fact that the bulk of 
the constraint is due to the LMC that has the largest 
velocity. The addition of the SMC either as independent 
tracer or as a second object in a pair to get the average 
center-of-mass values does not influence the constraint 
appreciably. 

The summary of all the mass constraints is presented 
in Table [1] In the left half of the table we show the con- 
straints obtained using host halos with any number of 
satellites with Vdrc > 50 km s^ 1 , while in the right half 
we show results obtained with the requirement that host 
halos host exactly two satellites with V C i rc > 50 km s _1 , 
as used in the study of Bll. The first column lists the 
constraints for the S sample, for the case where we mul- 
tiply the different likelihood distributions for each con- 
straint instead of computing the mass likelihood directly 
(S-multi and LGP-multi), the case when only properties 
of the LMC are used in the constraint, and the case where 
we assume that the SMC and LMC are a bound pair and 
correspond to a single subhalo with properties given by 
the mass-weighted average of the SMC and LMC prop- 
erties for the S and LGP samples (S-bound and LGP- 
bound). Columns 2 — 4 and 5 — 7 show the mass con- 
straints and the errors corresponding to the 68% and 
90% confidence intervals. M200 masses can be converted 
to M v i r using the conversion factor of M V i r /M 2 QQ — 1.21, 
computed directly the Bolshoi simulation for the halos of 
this mass range. 

4.3. Close satellite pairs 

Satellite pairs are quite rare: we find that only 2% of 
the MW-sized halos have a pair of subhalos with v c i rc > 
50km/s, separation < 40kpc, and relative velocity < 160 
km s _1 . This fraction drops to ~ 1/30000 for pairs more 
closely resembling the SMC-LMC pair. 

In this section we explore the incidence of the pop- 
ulation of close subhalo pairs and test whether their 
properties are distinct from the properties of the overall 
subhalo population. It is important to clarify these is- 
sues, because velocities and other properties of the SMC 
and LMC may be influenced by their mutual interac- 
tion (e.g.. iD'Onghia fc Lap 12001 iNichols et all 12011 
Bovlan-Kolchi n et al.l 120111 ). To select a sample of close 
pairs, we consider subhalos with v clrc > 50 km s within 
the MW-sized halos that have distance to the host cen- 
ter, mutual separation and absolute velocity similar to 
those of the MCs: 44 < r < 66 kpc, Ar < 60 kpc, 
150 < s < 400 km s -1 , and As < 300 km s _1 . We also 
estimate if the pairs are bound by computing the total 
energy assuming a two body system of two point masses 
with masses of the SMC and LMC as derived from the 
v c i rc — M relation in the simulation. Such approxima- 
tion should be viewed as a rough approximation for the 
fraction of bound systems in a given sample, but can be 
used as a reference. 

In figure [7] we show the separation and relative veloc- 
ities of the identified pairs (left panel) and velocity in 
the frame of the host halo of both pair members (right 
panel), with the ir-axis showing velocity of the pair mem- 
ber with the larger circular velocity. Red points indi- 
cate pairs that are bound according to our energy es- 



Gonzalez, Kravtsov & Gnedin 



TABLE 1 

Constraints on the halo mass of the Milky Way using different host samples and properties of the LMC and SMC 





Any number subs 


"circ > 50 km s 1 






Exactly two subs 


^circ > 50 km s 1 






Sample 


log(M 2O o/M ) 


68% c.i. 


90% 


c.i. 


log(Af 2 oo/M©) 


68% c.i. 


90% 


c.i 


S-sample 


12.06 


+0.08 -0.05 


+0.31 


-0.19 


12.03 


+0.06 -0.02 


+0.34 


-0.17 


S- multi 


12.24 


+0.13 -0.20 


+0.25 


-0.25 


12.15 


+0.15 -0.05 


+0.24 


-0.30 


LGP- multi 


12.44 


+0.07 -0.13 


+0.24 


-0.22 


12.26 


+0.14 -0.06 


+0.23 


-0.19 


LMC only 


12.18 


+0.10 -0.10 


+0.33 


-0.23 


12.17 


+0.10 -0.11 


+0.31 


-0.24 


S-bound 


12.13 


+0.13 -0.11 


+0.36 


-0.25 


12.09 


+0.11 -0.09 


+0.28 


-0.23 


LGP-bound 


12.17 


+0.16 -0.19 


+0.30 


-0.34 


11.99 


+0.16 -0.17 


+0.47 


-0.19 



Notes: columns 2 — 4 show the constraints for the case when host halos are allowed to have any number of satellites with u c i rc > 50 km 
s -1 , while columns 5 — 7 show the corresponding constraints for the case when host halos are restricted to have exactly two satellites with 
^circ > 50 km s — 1 . We present the mass errors corresponding to both the 68% and 90% confidence intervals from the derived likelihood 

distributions. 
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Fig. 7. — Properties of close subhalo pairs in MW mass halos. Left panel: pair separation and relative velocity; right panel: total speed 
relative to the host center for both pair members (horizontal axis for the member with larger circular velocity). Red points indicate subhalo 
pairs which are bound if we assume a two-body system. Blue rectangles show the MCs constraints (Ar < 60 kpc. As < 300 km s" 1 , 
150 < s < 400 km s _1 , 44 < ro < 66 kpc, and f c irc > 50 km s — 1 ). Green box in the left panel indicates very close pairs Ar < 15 kpc, 
As < 70 km s — 1 which contain a large fraction of fake pairs (see text for details). 



timate. The figure shows that only a small fraction of 
pairs has properties similar to those of the SMC-LMC 
pair Arjvic = 23.4 ± 10 kpc, and As M c = 128 ± 32 km 
s _1 (shown by the blue boxes in the two panels). The 
adopted range of separation of 23.4 ± 10 kpcis consider- 
ably larger than the actual observational error, which is 
of the same or der as the error in the ga lactocentric dis- 
tance: ~ 2 kpc (|Kallivavalil et al.ll2006al) . but it allows us 
to select MC-like pairs without imposing a prohibitively 
restrictive constraint on the subhalo pair configuration. 
Only « 60 satellite pairs out of 1140 outside the green 
region have separations and relative velocities similar to 
those of the SMC and LMC, ~ 90% of which are bound 
according to our criterion. 

The left panel of Fig. [7] shows a rather large number of 
pairs clustered at separation < 15kpc and As < 70 kpc. 
Most of these pairs appear to be artefacts due to fail- 
ures of the halo finder for subhalo candidates close to 



the resolution limit of the simulation. We have visu- 
ally inspected the dark matter density and velocity fields 
around a representative subset of these very close pairs 
and found that at separation 10 — 20kpc w 20% subhalo 
pairs do not have corresponding distinct peaks in dark 
matter density field or coherently moving clumps of dark 
matter (i.e., subhalos are fake). Some of the pairs ap- 
pear to be due to misidcntification, in which a subhalo 
undergoing a tidal disruption is identified as two subha- 
los with circular velocities close to the resolution limit. 
On the other hand, we find that for pairs with separation 
> 20kpc, almost all pairs have two clear distinct density 
peaks in the dark matter distribution and corresponding 
coherent velocity streams in the velocity field. 

If we add the constraints that the pair of subhalos in 
a host halos must have separation ArMC = 23.4 ± 10 
kpc, and Asmc = 128 ± 32 km s _1 , the resulting like- 
lihood for the Milky Way host halo mass is shown in 
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Fig. 8. — Likelihood distribution of the Milky Way halo mass 
using additional constraints that subhalo pairs have separation 
and velocity difference similar to those of the SMC-LMC pair: 
Ar M C = 23.4 ± 10 kpc, and As M C = 128 ±32 km s" 1 (black solid 
line). Likelihood distributions for subhalo pairs matching only the 
radial distance (orange line), circular velocity (red line), or rela- 
tive velocity with respect to the host halo (green line) are shown 
by the dotted lines. Less than 21% of close pairs are in halos of 
Af 2 oo > 1O 12 M . 

figure [51 Due to small number of close subhalo pairs sat- 
isfying the constraints, the mass likelihood distribution 
can be computed using only one of the MCs constraints 
at a time together with the separation and velocity dif- 
ference constraints. We find that < 21% of close pairs 
are in halos of Af 2 no > 1O 12 M The addition of the 
constraint on the pair separation and relative velocity 
thus pushes the MW halo mass constraint to considerably 
lower masses, as compared to the constraint with only the 
properties used by Bll. This illustrates that the actual 
constraint depends quite sensitively on which properties 
of the SMC-LMC system are chosen for the analysis, as 
is also clear from the large differences between mass like- 
lihood distributions for the individual properties shown 
in Fig. El 

5. DISCUSSION AND CONCLUSIONS 

In this study we have explored whether the environ- 
ment of the Milky Way affects properties of its two most 
massive satellites, the SMC and LMC. These satellites 
are rather rare around galaxies with the Milky Way lu- 
mi nosity and have rat her high velocities. As argued 
bv lBusha et~aTl (|2011al l. the properties of the SMC and 
LMC can be used for a useful independent estimate of the 
Milky Way virial mass. However, it would be straight- 
forward to interpret such estimate only if the Magellanic 
Clouds are not very unusual and their properties are not 
biased by the environment. 

To test for the effects of the environment, we compare 
velocity distributions for a sample of MC-sized subhalos 

6 This fraction is for the host halos restricted to have exactly 
two satellites with V c i rc > 50 km s — 1 . If we allow any number of 
satellites, this fraction increases to 26%. 



around Milky Way hosts selected simply by mass and a 
sample of such halos selected with additional restrictions 
on the distance to the nearest cluster and local galaxy 
density, designed to mimic the environment of the Local 
Group. We find that the velocity distribution of satellites 
in the latter sample is shifted somewhat to larger veloci- 
ties: in particular, the fraction of satellites with the LMC 
speed in the LGP sample is 23 ± 2% larger than in the 
S sample. Thus, the Local Group environment enhances 
the probability of high satellite velocities, although the 
effect is mild. 

We compute the likelihood distribution for the Milky 
Way halo mass usi ng the method similar to that of 
iBusha et al.l ((201 laT ) . and explore how this distribution 
depends on different properties of the SMC and LMC 
used as constraints. We find no significant effect of the 
environment on the mass estimate. We also find that the 
treatment of different MC properties as mutually inde- 
pendent does not bias the mass constraint. 

Using properties of the SMC an d LMC from the recent 
study bv iKallivavalil et al.l (120131 ). we derive constraint 
on the MW halo mass of log(M 20 o/M Q ) = 12.0^^ 
(90% confidence interval) for the S sample of halos se- 
lected without any environment restrictions. The mass 
constraint we derive is similar to that of Bll, even though 
the updated values of MC properties, such as velocities 
and their errors, arc quite different from the values used 
by Bll. The mass constraint is broadly consistent with 
other recent estimates of the Milky Way halo mass. 

The method of Bll does therefore appear to be a ro- 
bust way to measure the MW halo mass. A subtle issue, 
however, is that if some of the properties of the LMC or 
SMC used for the constraint are rare for the systems of 
the MW mass, the interpretation of the mass likelihood 
is not straightforward. The Milky Way in this case may 
be located in the tail of the distribution and interpre- 
tation of the peak of the likelihood as most likely mass 
of the MW halo is not correct. It remains to be seen 
how rare particular properties of the LMC and SMC are. 
We do know, for example, that the probability of having 
two satellites of the SMC and LMC luminosity is by it- 
self quite rare in the hosts of the MW luminosity (e.g., 
ILiu et al.ll2^TTbl:IRobotham et al.|[20Tl . 

LMC and SMC are not just two unrelated luminous 
satellites: they are a close pair. We show in this study 
that close satellite pairs are quite rare: pairs with masses 
and separations similar to those of the LMC and SMC 
are found only in one out of w 30000 MW-sized halos. 
Obse rvations also indicate that such close pairs are very 
rare ([ James fc Ivorvl [20111) . We find that satellites in 
most close pairs with properties similar to the MCs are 
likely to be bound to each other. Interestingly, the halo 
mass likelihood distribution for host halos constrained to 
have MC-like close pairs of subhalos (Ar^ic = 23.4 ± 10 
kpc, and Asmc = 128 ±32 km s _1 ) is quite different from 
the global likelihood from which the MW halo mass con- 
straint discussed above was derived. In particular, less 
than 21% of host halos with M 20 o > 10 12 M Q host MC- 
like close pairs. Taking into account the close separation 
of the MCs in the Bll method results in the shift of the 
MW halo mass estimate to smaller masses with the peak 
shifting approximately by a factor of two (see Fig. |SJ. 
The reason for this shift is the fact that in smaller ha- 
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los it is more likely to get a pair with small separation 
and relatively small velocity difference by chance (not 
necessarily bound) , mainly because velocity dispersion is 
smaller in smaller mass halos. 

This example clearly shows that a great care should 
be taken in choosing which of the satellite properties are 
used for the MW mass constraint. 
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